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ABSTRACT 


A  local  deliver?  aodel  was  developed  for  a  repair  facil¬ 
ity-stock  point  system,  given  one  or  aore  supported  produc¬ 
tion  lines  and  each  coaponent  repaired  nay  require  aore  than 
one  part.  Both  deterainistic  and  randoa  deaands  were  con¬ 
sidered.  The  objective  function  was  total  expected  trans¬ 
portation  and  delay  costs  per  day.  In  the  deterainistic 
case  the  total  cost  curve  was  discontinuous  and  the  optiaal 
delivery  policy  could  only  be  determined  by  exhaustive  enu¬ 
meration.  A  computer  simulation  model  was  needed  for  the 
random  demand  case.  The  siaulation  aodel  was  also  extended 
to  allow  randoa  issue  processing  tine  and  a  remote  warehouse 
sited  close  to  the  repair  facility.  The  results  of  the 
siaulation  showed  that  point  of  entry  effectiveness  and 
non-local  response  tines  were  key  factors  of  expected  delay 
costs  and  that  these  costs  could  be  reduced  through  the  use 
of  a  reaote  warehouse,  aore  importantly,  providing  the  best 
support  to  custoaers  requiring  the  fewest  parts  per  conpo- 
nent  repaired  will  give  the  ainiaua  expected  delay  cost. 
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I.  I WTRODOCTION 


If  all  required  materials  were  available  at  the  right 
time  and  place  for  a  reasonable  price,  no  manager,  business, 
or  government  agency  would  choose  to  stock  them.  Unfortu¬ 
nately  this  is  not  the  case  and  both  the  Department  of 
Defense  and  the  Navy  maintain  large  stocks  of  material  in 
support  of  their  missions.  With  increasingly  complex  and 
specialized  weapon  systems,  the  sources  of  supply  are  becom¬ 
ing  more  scarce  and  procurement  lead  times  are  increasing, 
resulting  in  the  need  for  increased  range  and  depth  of  sup¬ 
port,  Meanwhile  pressures  to  decrease  the  federal  budget 
deficit  and  a  high  inflation  rate  have  often  forced  the 
Navy,  as  well  as  other  government  agencies,  to  operate  on 
budget  allotments  which  may  be  declining  in  purchasing 
power.  To  maintain  previous  levels  of  service,  increases  in 
operational  efficiency  and  worker  productivity  at  least 
equal  to  that  being  obtained  by  private  industry  are 
required. 

The  consolidation  of  support  facilities  within  the  Navy 
has  been  one  method  of  improving  efficiency.  The  develop¬ 
ment  of  centralized  Inventory  Control  points  (ICPs)  have 
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certainly  had  significant  impact  on  the  supply  system. 
Through  the  collection  and  manipulation  of  a  system-wide 
data  base,  more  intelligent  provisioning,  outfitting,  budg¬ 
eting,  and  stocks ge  decisions  have  been  possible.  Providing 
vorld-vide  asset  visibility  and  centralized  procurements 
have  also  offered  improved  support  at  a  reduced  cost.  It  is 
expected  further  improvements  will  still  be  made  in  this 
area  in  the  future  [Hef.  1], 

Such  of  the  success  of  the  ICP  effort,  however,  has  to 
be  attributed  to  the  development  of  high  speed  communica¬ 
tions  systems  used  to  transfer  information  to  the  ICP  and 
the  development  of  high  speed  and  high  capacity  computers 
and  peripherals  to  process  this  information,  without  the 
necessary  capital  investment  in  the  above  productivity 
enhancing  systems,  the  ICP  would  likely  be  a  much  less 
effective  and  desirable  entity. 

Consolidations  have  been  occurring  in  other  areas  as 
well.  flajor  stock  points  at  Newport,  Rhode  Island  and  Long 
Beach,  California  have  essentially  been  closed  or  consoli¬ 
dated  with  other  support  activities.  Material  for  fleet 
issue  has  been  consolidated  at  regional  Naval  Supply  Centers 
(NSCs)  located  at  major  demand  sources.  The  most  recent 
moves  have  been  to  consolidate  wholesale  supply  support  for 
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several  Naval  Air  Rework  Facilities  (NARFs)  at  nearby  Naval 
Supply  Centers.  Previous  support  had  been  provided  by  Naval 
Air  Station  supply  departments  where  those  NARFs  are 
located.  Since  the  supply  centers  often  carry  material 
under  Defense  Logistic  Agency  (DLA)  funding  as  well  as  that 
provided  by  the  Navy  Stock  Fund  and  Navy  Industrial  Fund, 
stock  range  and  depth  should  improve  over  that  which  was 
previously  available  at  the  air  station.  This  improved 
stock  position  should  lead  to  improved  point  of  entry  (POE) 
effectiveness  and  thus  improved  customer  support,  other 
things  being  equal.  These  consolidations  of  support  are 
made  economically  more  attractive  when  the  supply  centers 
install  capital  intensive,  productivity  enhancing  automated 
material  handling  systems  such  as  NISTARS  (Naval  Integrated 
Storage  and  Retrieval  System). 

However,  by  centralizing  material  at  regional  centers, 
distances  that  material  must  move  after  issue  to  reach  the 
customer  may  increase  substantially.  Not  only  would  this 
possibly  increase  transportation  costs,  but  more  importantly 
it  would  likely  delay  the  receipt  of  required  parts  on  the 
customer's  production  line,  With  components  under  repair 
awaiting  parts,  either  test  bench  or  shop  space  is  occupied 
or  maintenance  time  must  be  used  instead  to  consolidate  the 
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pieces  of  the  component  in  progress  and  store  them  together 
until  the  required  parts  are  received.  In  either  case 
valuable  production  resources  are  lost,  thus  incurring  some 
delay  cost. 


II .  PRSVioos  §122125 


The  first  Naval  Air  Rework  Facility  (NARF)  wholesale 
support  consolidation  was  that  of  NARP  Alameda  and  Naval 
Supply  Center  (NSC)  Oakland,  which  occurred  in  October  1979. 
Prior  to  that  consolidation.  Grant  [Ref.  2]  attempted  to 
quantify  the  production  delay  costs  caused  by  not  having 
repair  parts  immediately  obtainable  when  needed  by  research¬ 
ing  NARF  Alameda  procedures  and  records.  Although  some 
costs,  such  as  cannibalizat ions ,  had  avenues  for  documenta¬ 
tion,  others  did  not  and  he  was  unable  to  develop  a  firm 
relationship  between  delivery  times  and  delay  costs.  In  the 
preparation  of  his  thesis  Grant  conducted  interviews  at  NARF 
Alameda  and  many  shared  one  common  view.  Overall  availabil¬ 
ity  is  much  more  important  than  the  rapid  delivery  of  less 
than  all  the  parts  required  to  repair  a  component. 

It  is  really  the  slowest  delivery  which  sets  the  pace  of 
the  repair  action  and  should  be  used  to  determine  production 
delay  costs  caused  by  the  lack  of  repair  parts.  For  exam¬ 
ple,  if  four  parts  were  required  and  three  were  delivered  in 
twenty  minutes  and  the  fourth  was  not  delivered  until  two 
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weeks  later,  the  component  (barring  cannibalization)  would 
spend  two  weeks  awaiting  repair  parts. 

In  a  second  thesis  Davidson  [Ref.  3]  conducted  an  analy¬ 
sis  of  three  direct  delivery  models  which  were  initially 
proposed  by  acHasters  [Ref.  4].  These  models  were  based  on 
a  single  customer  (such  as  the  jet  engine  repair  line  at 
WAR?  Alameda)  and  a  single  repair  part  which  may  need  to  be 
replaced  and  thus  require  requisitioning  for  each  inducted 
component.  The  demand  for  this  part  was  considered  a  Ber¬ 
noulli  trial  with  a  fired  probability  of  demand  (p)  for  each 
induction. 

HcMasters  and  Davidson  attempted  to  minimize  expected 
costs  where  total  cost  was  the  sum  of  transportation  cost  (a 
fixed  charge  per  delivery)  and  delay  costs  (a  fixed  charge 
per  component  per  unit  of  time  delayed  due  to  the  lack  of 
the  repair  part) .  The  only  delays  considered  were  those 
caused  by  the  transportation  system  (i.e.  material  avail¬ 
ability  was  not  considered)  and  the  unit  of  time  was  defined 
as  the  time  between  component  inductions  on  the  repair  line. 
Expected  total  costs  were  calculated,  but  due  to  analytic 
complexities  of  these  models,  closed  form  optimizations  for 
the  models  were  not  possible.  Instead,  a  parametric  analy¬ 
sis  was  conducted  for  each  of  the  three  delivery  plans. 
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Davidson  shoved  that,  although  ths  plans  considered  differed 
significantly  in  form  and  emphasis,  there  was  little 
difference  in  the  optiaal  expected  costs  for  each.  She  also 
shoved  that  varying  the  delay  cost  per  period  (CD)  had  a 
auch  greater  iapact  on  the  optiaal  total  cost  than  varying 
the  paraaeter  p,  although  increasing  p  did  increase  cost. 

This  thesis  will  extend  the  work  of  References  3  and  4. 
It  will  consider  systeas  with  one  or  aore  customers,  each  no 
longer  limited  to  one  repair  part  per  induction.  Chapter  3 
broadly  suaaarizes  the  earlier  models  and  then  discusses 
additional  assumptions  needed  to  generalize  these  models. 
Finally  it  presents  a  new  model  for  the  deterministic  case. 
Chapter  4  studies  stochastic  versions  of  the  new  model  and 
Chapter  5  considers  the  impacts  on  this  model  of  locating 
aaterial  at  the  customer's  site.  Chapter  6  presents  a  sum¬ 
mary  and  conclusions. 


III.  SUMMARY  OF  MODELS.  CONSIDERED 


The  basic  system  being  modeled  is  diagrammed  in  Figure  1 


Failed  Repairables 


Requisitions 


\ 


\ 


|Repair  Facility) 

/  \ 

Repaired  Components  Repair  Parts 


Stock  Point 


Pigure  1:  Customer-Stock  Point  Relationship 


above.  The  industrial  customers  considered,  such  as  produc¬ 
tion  lines  at  a  Naval  Air  Rework  Facility,  induct  components 
for  repair,  troubleshoot  each  component,  requisition  any 
required  repair  pacts,  and,  upon  receipt  of  those  parts, 
complete  repair  of  the  failed  item.  Earlier  studies  [Ref.  3 
and  U]  considered  alternative  transportation  systems  for 
delivering  a  given  required  repair  part  from  the  stock  point 
to  the  customer  and  attempted  to  minimize  the  sun  of 
expected  transportation  costs  and  expected  customer  delay 
costs.  McMasters  [Ref.  4]  also  addressed  the  establishment 
of  an  On-Site  Inventory  System  (OSIS)  at  the  customer’s 
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location  to  expedite  delivery  and  reduce  customer  delay 
costs.  This  study  will  address  the  OSIS  in  a  later  chapter. 

McMasters  proposed  three  basic  local  delivery  options. 
These  were: 

t.  Deliveries  are  made  at  the  end  of  N  periods  if  there 
has  been  at  least  one  demand  during  that  time  frame. 

2.  Deliveries  are  made  as  soon  as  K  issues  accumulate. 

3.  A  delivery  is  made  in  the  (N-I)st  period  after  the 
first  demand  following  a  delivery. 

Initially  this  study  will  consider  only  Option  1. 

Davidson  [ Hef .  3]  shows  that  for  the  single  customer  case 
all  3  models  display  nearly  equal  cost  structures  and  recom¬ 
mended  Option  1  as  a  quite  reasonable  strategy.  Option  1 
also  seems  best  suited  to  non-industrial  activities  (such  as 
ships  in  port)  who  must  schedule  workers  based  on  parts 
availability.  By  knowing  when  deliveries  are  made,  requisi¬ 
tion  status#  and  the  ship's  operating  schedule,  supervisors 
can  estimate  when  technicians  must  work  extra  hours  or  when 
they  can  be  given  extra  time  off.  Also#  by  knowing  the 
delivery  schedule#  extraordinary  action  can  be  taken  if  sys¬ 
tem  response  will  not  satisfy  a  particularly  critical  need. 
Option  1  is  also  representative  of  the  way  many  stock  points 
currently  operate  their  local  delivery  system. 
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This  paper  will  modify  the  above  aodel  so  that  costs  for 
a  multi-customer  (or  multi-production  line) ,  aulti-item 
inventory  system  can  be  considered,  cost  structures  of  the 
aodel  will  be  studied  in  the  hope  of  deteraining  rules  for  a 
cost  minimizing  delivery  plan.  The  effect  cf  locating  sup¬ 
porting  stocks  at  the  site  of  the  customer  will  also  be 
studied. 

A.  ASSUMPTIONS  MADE  IN  GENERALIZING  FROM  PREVIOUS  MODELS 

In  generalizing  to  multi-item,  multi-customer  systems, 
some  assumptions  from  the  earlier  aodel  aust  be  modified  and 
some  additional  assumptions  made.  First,  McHasters*  aodel 
referred  to  a  time  period  as  "the  time  between  component 
inductions  on  the  production  line".  This  is  not  convenient 
for  the  multi-production  line  environment  where  different 
customers  may  have  different  periods  between  inductions. 
Therefore  a  common  denominator  for  time  among  all  production 
lines  or  customers,  the  work  day,  is  used  in  this  study  as 
the  unit  for  time. 

As  with  the  earlier  aodel,  transportation  costs  will  be 
considered  as  a  fixed  charge  per  shipment.  In  a  multi-item 
inventory  it  might  be  more  realistic  to  allocate  charges  by 
weight  or  volume,  particularly  if  deliveries  were  con¬ 
strained  by  one  of  those  parameters.  This  was  not  done. 
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however,  since  local  custoaers  are  being  considered  and  it 
is  felt  that  local  deliveries  are  not  usually  capacity 
constrained.  Moreover,  by  not  being  forced  to  specify 
specific  itea  weight  and  cube,  the  model  could  reaain  more 
generally  applicable. 

Next  it  will  be  assuaed  that  all  requisitions  are  homo¬ 
geneous  within  the  issue  and  transportation  system.  This 

aeans  requisitions  are  distinguished  by  requisition  number 

/ 

and  customer  only,  and  not  by  priority,  weight  or  cube,  or 
itea  required.  Although  requisition  quantity  may  be  greater 
than  one,  issue  of  partial  quantities  is  not  considered. 
Finally,  all  requisitions  are  assuaed  to  receive  the  same 
processing  within  the  system. 

By  making  these  assumptions  the  multi-repair  part  local 
delivery  problem  becomes  one  of  tracking  multiple  requisi¬ 
tions  for  each  itea  under  repair.  Although  this  simplifica¬ 
tion  does  not  allow  for  interdependent  response  times,  such 
as  aight  be  expected  when  spares  are  driven  to  a  not-in- 
stock  position,  it  does  allow  for  interdependence  in  the 
transportation  system  for  the  "ship  every  K  issues"  case. 

It  will  be  assuaed  that  component  inductions  are  made  a 
fixed  period  apart.  This  period  is  a  given  parameter  7  for 
each  custoaer.  Although  it  is  usually  determined  by  the 


nuaber  of  components  scheduled  for  repair  in  the  current 
calendar  quarter,  it  can  also  be  considered  the  aaintenance 
time  required  for  component  repair  given  a  aaintenance 
resource  allocation.  As  more  components  are  required  per 
quarter,  Y  will  decrease  and  the  shop  supervisor  will  have 
to  assign  more  production  resources. 

In  the  proposed  model  it  is  assumed  the  actual  repair 
period  is  divided  into  three  major  phases.  The  first  is  the 
troubleshooting  phase  which  is  hypothesized  to  take  one  half 
the  scheduled  maintenance  time,  or  Y/2  days.  During  this 
time  the  component  is  disassembled  and  all  the  parts  which 
need  to  be  replaced  are  determined. 

Phase  two  of  the  repair  process  is  the  "obtain  the 
repair  parts"  phase.  It  consists  of  ordering  all  required 
parts  and  waiting  for  their  receipt.  Since  all  requirements 
were  determined  in  phase  one,  this  phase  takes  essentially 
no  maintenance  effort.  During  this  time,  aaintenance 
resources  can  be  allocated  to  other  jobs.  Delay  in  receiv¬ 
ing  the  required  parts  does  incur  costs  in  work  in  process 
inventory,  maintenance  test  bench  space  occupied,  and  inef¬ 
ficiencies  cause  by  moving  maintenance  personnel  between 
jobs.  For  that  reason,  delay  charges  are  assumed  to  be 


assessed  at  a  fixed  rate  (CD  dollars)  per  component  per  day 
for  the  tioe  spent  awaiting  repair  parts. 

The  final  phase  of  the  repair  process  is  the  actual 
repair  of  the  component.  This  includes  replacing  the  failed 
parts,  component  reassembly,  and  final  test.  This  last 
phase  is  allocated  50*  of  the  maintenance  effort,  or  Y/2 
days 

It  may  be  somewhat  confusing  that  Y  does  not  equal  the 
sum  of  the  three  phases  of  the  repair  cycle.  This  is 
because  more  than  one  component  can  be  in  process  at  any  one 
time,  and  should  be  if  a  component  is  awaiting  repair  parts. 
Y  is  the  time  between  inductions,  the  average  time  between 
repair  completions,  and,  in  this  model,  the  time  between 
submission  of  requisition  batches.  The  average  component 
turn  around  time  is  the  sum  of  the  tine  spent  in  each  of  the 
repair  phases,  or  Y  plus  average  delay  time. 

The  above  three-phase  repair  process  assumes  that  all 
the  parts  required  can  be  determined  in  phase  one  and 
ignores  any  parts  broken  or  discovered  defective  during 
reassembly.  This  is  considered  realistic  since  the  rework 
facility  has  typically  been  repairing  the  component  in  ques¬ 
tion  for  a  long  tine  and  these  last  minute  demands  can  often 
be  anticipated. 
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By  assessing  delay  cost  at  a  constant  rate  (CD)  until 
all  parts  are  received,  any  benefits  of  receiving  some  but 
not  all  repair  parts  are  ignored.  This  does  suppress  any 
benefits  from  cannibalization,  but  cannibalization  costs  can 
be  high  and  the  above  does  seen  the  most  fair  way  to  levy 
delay  costs. 

B.  DERIVING  THE  EXPECTED  COST  FUNCTION  FOR  THE  IMMEDIATE 

ISSUE  CASE 

Consider  the  single  customer  case  where  the  system  has  a 
transportation  cost  of  CT  dollars  per  delivery,  and  the  cus¬ 
tomer  has  a  delay  cost  rate  of  CD  dollars  per  day  per  compo¬ 
nent  and  a  scheduled  induction  period  of  I  days.  The 
decision  variable  for  the  system  is  N,  the  periodicity,  in 
days,  of  deliveries.  The  objective  will  be  to  minimize  the 
average  daily  total  cost  where 

Average  Total  =  Average  Transportation  ♦  Average  Delay 
Daily  Cost  Cost  Per  Day  Cost  Per  Day 

or 

ADC  (N)  ■  TC  ♦  DC 

To  derive  the  average  total  daily  cost,  the  process  must 
be  examined  a  little  nora  closely.  Consider  first  the  sin¬ 
gle  customer  case.  As  long  as  N,  the  cumber  of  days  between 
deliveries,  and  7,  the  days  between  inductions  of  a  compo¬ 
nent  for  repair,  are  rational,  this  will  be  a  renewal 
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process.  If  rational,  uY*vN  for  some  integers  u  and  v  and 
the  system  will  cycle  every  u  inductions  or  v  deliveries. 

To  determine  long  run  time-average  costs,  costs  will  only 
have  to  be  averaged  over  a  cycle.  In  the  case  of  deliveries 
which  cost  CT  dollars  each,  the  total  cost  for  the  v 
deliveries  of  the' cycle  would  be  vCT.  Since  deliveries  are 
M  days  apart,  the  total  length  of  the  cycle  is  vN  days. 
Dividing  the  total  delivery  cost  per  cycle  by  the  days  per 
cycle,  average  daily  transportation  cost  becomes 
TC  »  _v_  *  C|T 

Delay  costs  are  a  little  bit  more  complex  for  they  are  a 
function  of  both  H  and  T.  In  all,  three  different  parameter 
conditions  can  be  considered.  First  consider  delay  costs 
when  H  is  less  than  T.  This  implies  deliveries  are  more 
frequent  than  inductions  on  the  one  production  line  consid¬ 
ered.  Although  this  may  seem  unrealistic  in  the  single 
customer  case  since  some  deliveries  would  consist  of  no 
requisitions,  it  could  easily  arise  when  multiple  customers 
at  a  single  location  or  on  a  single  local  delivery  route  are 
considered.  In  any  case,  Pigure  2a  shows  the  time  until  the 
next  delivery  for  a  delivery  schedule  with  N  equal  to  4 
days.  Superimposed  on  the  x-axis  and  marked  with  triangles 
are  the  times  when  the  requisition  submissions  would  take 
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place  if  the  induction  periodicity,  7,  equalled  4.5  days  and 
the  first  delivery  and  order  were  concurrent,  As  can  be 
seen  from  the  figure,  the  delay  for  the  first  and  ninth 
inductions  would  be  the  sane  and  thus  as  long  as  N  and  T 
renain  constant,  the  length  of  component  delay  would  cycle 
every  eight  inductions.  Shown  in  the  bottom  graph  of  Pigure 
2  (Pigure  2b)  is  the  delay  in  days  for  each  component.  Note 
that  if  the  initial  delivery  were  a  bit  later  it  would 
increase  the  delay  time  for  each  of  the  seven  subsequent 
induction  in  the  cycle.  Thus  when  calculating  average  com¬ 
ponent  delay,  this  phase  factor,  call  it  q,  based  on  initial 
conditions,  should  be  added.  However,  it  should  be  obvious 
that  any  optimal  delivery  plan  should  have  initial  condi¬ 
tions  adjusted  so  that  this  q  would  be  equal  to  zero.  For 
this  reason  q  will  be  assumed  zero  for  the  rest  of  this 
study. 

Appendix  A  assumes  both  N  a nd  I  rational  and  solves  for 
the  values  of  u  and  v  mentioned  in  the  above  renewal  process 
argument.  It  derives  component  delay  as  a  recursion  rela¬ 
tion  and  shows  that  the  number  of  inductions  in  the  cycle  is 
N/L,  where  L  is  the  largest  real  number  common  to  both  N  and 
I.  L  is  defined  such  that  Y/L  and  N/L  are  both  integers, 
integers  which  are  actually  the  u  and  v  which  were  referred 
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to  earlier.  The  appendix  then  goes  on  to  show  that 
component  delay  accepts  N/L  evenly  stepped  discrete  values 
and  the  mean  of  these  values  is 

Average  Delay  per  Component  =  ( (N  -  L)/  2)  , 

yielding  an  average  delay  cost  of 

Average  Delay  Cost  per  Component  =  CD  ( (N  -  L)/  2)  . 
Since  this  is  the  average  delay  cost  per  component  and  u 
components . were  inducted  per  cycle,  the  total  delay  cost  per 
cycle  would  be  uCD(N-L)/2.  The  length  of  the  cycle  would  be 
ul  days  so  the  average  daily  delay  cost  becomes 

DC  =  qi  fN-H  (3.1) 

Next  consider  the  case  where  N=T,  or  where  deliveries 
and  inductions  have  the  same  periodicity.  In  this  case  all 
components  would  experience  the  same  delay.  As  long  as 
deliveries  and  orders  were  perfectly  phased  (i.e.  initial 
conditions  were  right),  each  component  would  experience  zero 
delay  and  hence  zero  delay  cost.  Note  that  this  zero  delay 
cost  would  be  provided  by  equation  (3.1)  since  L  would  be 
equal  to  N. 

Lastly,  consider  the  case  where  N  is  greater  than  I.  A 
special  case  of  this  condition  is  N=iT  for  some  integer  i 
greater  than  1.  Under  this  condition  all  the  deliveries 
will  still  be  at  the  same  point  in  each  repair  cycle  but  now 
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■ore  than  one  component  will  be  awaiting  repair  parts. 
Assuming  cost  minimizing  initial  conditions,  one  component 
would  experience  no  delay.  Since  i  components  would  have 
been  inducted  since  the  last  delivery,  (i-1)  components  must 
have  been  waiting  repair  parts  the  last  Y-day  induction  per¬ 
iod,  (i-2)  the  induction  period  before  that,  and  so  on. 

Thus  the  total  component  delay  per  shipment  would  be 
I((i-1)  +  {i-2)+...  +  H0),  or 
*(i(i^D) 

Since  there  were  v  shipments  per  renewal  cycle  and  N  days 
between  shipments,  the  average  daily  delay  cost  becomes 
DC  =  _j_-  v-CD*  Yi  (i^l)  -  CJ, 

But  H*iY  so 

DC  »  CD  fi-11  (3.2) 

This  is  the  same  deterministic  delay  cost  eguation  as  was 
developed  by  HcHasters.  note  that  if  i=1,  delay  costs  are 
zero  as  was  predicted  earlier. 

Eguation  (3.2)  is  also  a  degenerate  form  of  equation 
(3.1)  for  the  special  case  N=iY.  Since  N=iY,  L  must  assume 
the  value  Y  as  long  as  i  is  integer.  Osing  this  fact,  egua¬ 
tion  (3.1)  becomes 
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DC  * 


*  CD  (j  -  .1) 


Next  consider  the  general  case  where  N  is  greater  than 
Y,  or  deliveries  are  less  frequent  than  inductions.  Alt¬ 
hough  there  may  be  more  than  one  component  awaiting  repair 
parts  at  any  one  time,  steady  state  average  daily  delay 
costs  can  still  be  obtained.  Approaching  the  problem  in  a 
method  similar  to  the  N  less  than  Y  case.  Figure  3a  is  a 
graph  of  the  delivery  schedule  for  N  equals  5  and  Y  equals  3 
days.  The  requisition  times  for  a  customer  are  marked  as 
triangles  on  the  abscissa.  Note  that  in  the  case  illus¬ 
trated  delay  times  within  the  cycle  are  not  monotone 
decreasing  as  they  were  in  the  previous  case  (Figure  2)  ,  but 
delays  are  still  in  multiples  of  L.  As  the  derivations  in 
Appendix  A  still  hold,  average  delay  cost  reduces  to  equa¬ 
tion  (3.1)  again,  or 

°C  =  ^  |L5_l) 

Combining  the  transportation  cost  and  delay  cost  terms, 
the  overall  single  customer  daily  cost  function  becomes 

(*-H  .  <3-31 

where  L  is  the  largest  number  such  that  Y/L  and  N/L  are 
integers . 


ADC  (N) 


CJ.  *  ££ 
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For  the  two  customers  case  where  custoner  parameters  are 
CD1,  r  1,  01,  and  CD2,  Y2,  D2  respectively,  delay  costs  for 
the  customers  are  summed  to  obtain  the  system  delay  cost  or 
ADC  m  a  eg  »  Cjl  .  |H  -  Li  j  ♦  Cttj.jN  -  L2|  .  (3.4) 

Equation  (3.4)  does  assume  that  N,  II,  and  Y2  are  rational, 
or  that  deliveries  can  be  phased  such  that  at  one  point  in 
time  both  customers  can  experience  zero  delay. 

By  ignoring  the  L  terms,  an  upper  bound  approximation 
can  be  obtained  for  (3.4)  .  Generalizing  this  approximation 
to  n  customers  the  average  total  daily  cost  can  be  expressed 
as 

ADC(N)  =  CT  ♦  ♦  ...  ♦  Cf&j  (3.5) 

C.  OPTIMIZING  THE  AVERAGE  COST  FUNCTION 

Even  though  the  cost  expression  given  by  equation  (3.3) 
is  only  for  the  deterministic  case,  it  is  not  easily  minim¬ 
ized.  The  term  which  cause  the  difficulty  in  optimization 
involves  L,  which  is  not  continuous  in  N.  Kith  that  being 
the  case,  one  way  to  "optimize"  the  function  is  to  compute 
costs  for  the  various  values  of  N  which  are  of  interest  and 
select,  as  optimal,  that  N  which  gives  minimum  plotted  cost. 
Before  doing  this,  however,  it  is  possible  to  get  an  upper 
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bound  on  costs  by  deleting  the  L  term  froa  the  cost  equa¬ 
tion.  The  resulting  approx  iaation  is  continuous  in  N  and 
can  be  ainiaized  using  the  calculus.  Using  the  aore  general 
n  custoaer  case  or  equation  (3.5)  and  solving  for  the  first 
order  conditions  for  minimization. 

Si22' IP  **•(=?**=?**•••  *=»)■• 
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CD1  ♦  CD2  ♦  ... ♦  CDn 
II  12  In 


(3.6) 


Checking  the  second  order  conditions 

m*  -  f* 


which  is  greater  than  0  for  positive  N  and  CT,  and  thus  N* 
given  by  equation  (3.6)  ainiaizes  (3.5). 

Figure  4  investigates  the  shape  of  this  bounding  cost 
function  for  the  single  customer  case  with  CD=100  dollars, 
CT*100  dollars,  and  1*3  days.  It  shows  the  total  average 
cost  and  its  components,  transportation  cost  and  delay  cost 
plotted  for  various  N,  the  delivery  periodicity.  When  n*1, 
equation  (3.5)  is  similar  to  the  Hadley  and  Whitin  [Ref.  5] 
Deterministic  Lot  Inventory  Hodel  cost  function.  Ks  a 
consequence,  the  square-root  formula  for  N*  resembles  that 
of  the  economic  order  quantity. 
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DELIVERY  PERIODICITY  (DAYS) 


Figure  4:  Approximate  Cost  Punction  components 

Figure  5  offers  a  comparison  of  the  upper  bound  function 
and  the  exact  cost  functions  in  both  the  one  and  two  cus¬ 
tomer  cases.  The  top  graph  is  that  for  a  customer  hawing 
CD*100  and  T*7  and  a  system  CD*10Q.  While  the  upper  bound 
(approximate)  curve  is  smooth  and  has  a  minimum  near  3.74, 
the  exact  costs  as  derived  from  equation  (3.3)  and  plotted 
as  triangles,  would  not  have  a  smooth  curve.  Although 
transportation  costs  are  decreasing  as  N  is  increased,  the 


31 


3  INGLE  CUSTOMER 


delay  cost  tera  decreases  and  increases  depending  on  rela¬ 
tive  values  of  L.  At  N=7  delay  costs  go  to  zero,  giving  a 
ainiaua  total  cost.  Note  however  that  varying  N  a  small 
amount  fro  a  this  point  gives  costs  which  are  near  the  upper 
bound  function  since  L  decreases  sharply.  Since  the  cost 
function  is  not  continuous,  lines  should  not  be  drawn  bet¬ 
ween  these  exact  cost  points.  Such  lines  would  only  encour¬ 
age  interpolation  which  could  lead  to  invalid  conclusions. 

The  lower  graph  in  Figure  5  is  for  the  two  customers 
case.  The  second  customer  is  assumed  to  have  CD=3 0  and  Y*3. 
Again  the  exact  cost  points,  this  time  from  equation  (3.4), 
are  plotted  as  triangles  and  again  the  ainiaua  cost  is  not 
necessarily  near  the  ainiaua  of  the  approxiaate  function. 
This  tine  the  exact  cost  points  are  the  sun  of  three  terms 
in  the  cost  equation  which  act  seemingly  independently. 
Although  transportation  costs  are  decreasing  monotonies lly 
as  N  is  increased,  the  two  delay  cost  terns  increase  and 
decrease  depending  on  the  values  for  LI  and  L2.  Note  that 
at  N*3  the  delay  cost  for  Customer  2  goes  to  zero  since  N*r, 
but  the  delay  costs  for  Customer  1  get  a  much  smaller  break. 
LI  at  H*3  is  1  so  Custoaer  1  delay  costs  are  2/3  of  the 
upper  bound  amount.  Also  note  that  very  slight  changes  froa 
N*3  (such  as  N*3.001)  give  very  snail  values  for  Li  and  L2 
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and  thus  delay  costs  are  very  near  the  upper  bound  amount. 
Again  the  apparent  minimum  in  Figure  5  appears  at  N«7.  At 
this  point  Customer  1  experiences  zero  delay  costs  and  Cus¬ 
tomer  2's  delay  cost  is  6/7  of  the  worst  case  amount.  This 
would  not  be  the  minimum  except  for  the  fact  that  the  delay 
cost  rate  for  Customer  1  (100)  is  significantly  greater  than 
that  of  Customer  2  (30) .  Although  not  plotted,  again  minor 
variations  either  side  of  S»7  yield  delay  costs  and  total 
costs  near  the  worst  case  curve. 

If  N  is  greater  than  7,  S  is  greater  than  T  for  both 
customers  and  both  customers  will  experience  some  delay 
costs.  Even  though  transportation  costs  are  decreasing,  it 
appears  this  decrease  is  less  that  the  increase  in  delay 
costs  and  H»7  is  the  true  minimum.  Also  note  that  in  the 
two  customers  case  the  exact  cost  points  more  closely 
approximate  the  upper  bound  case  in  shape  (although  this 
approximation  is  still  quite  poor).  As  more  customers  are 
added  to  the  system,  more  delay  costs  are  added  to  the  total 
cost  expression.  Thus  each  individual  delay  cost  term  is  a 
smaller  proportion  of  total  costs  and  as  long  as  the  I 
values  are  not  the  sane,  the  upper  bound  approximation 
should  improve. 
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Fortunately,  values  for  N  are  not  chosen  in  a  continuous 
■anner  and  deliveries  are  usually  aade  every  half  day,  day, 
or  soaething  like  that.  If  Y  and  CO  values  for  all  custoa- 
ers  are  known  exactly,  the  total  average  cost  for  each  value 
of  S  can  be  calculated  and  that  generating  the  minimum  costs 
would  be  chosen  as  optiaal.  If  only  approximate  values  for 
t  are  known,  perhaps  using  the  approximate  or  worst  case 
function  would  be  the  best  strategy. 
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IV .  SIMULATION  MODELS 


As  can  be  seen  froi  the  previous  chapter,  deriving 
delivery  plans  that  minimize  cost  is  difficult  at  best  even 
when  using  a  relatively  siaple  deterministic  model.  when 
coaplicating  factors  such  as  stochastic  demands,  lead  times, 
or  induction  periods  are  included,  the  mathematics  quickly 

becomes  extremely  complex  and  is  not  easily  analyzed  through 

| 

the  use  of  the  calculus.  Por  this  reason  a  simulation  model 
of  the  system  was  written  in  the  Simscript  II. 5  language. 
This  program  is  an  event  step  simulation  and  a  listing  of 
the  basic  program  is  included  as  Appendix  B. 

In  an  event  step  simulation,  specific  events  are  sche¬ 
duled  and  executed  at  specific  points  in  time.  These  events 
often  lead  to  other  events,  which  are  then  scheduled  during 
execution.  Figure  6  is  a  broad  flowchart  of  the  main  events 
used  to  determine  cost  estimates  for  the  system  under  study. 
The  simulation  allows  using  either  the  ship  every  N  days  or 
the  ship  every  K  requisitions  delivery  options.  If  the  ship 
every  K  requisitions  option  is  used,  each  time  an  issue  is 
made  the  program  determines  if  K  requisitions  have  accumu¬ 
lated.  If  so,  a  delivery  is  scheduled.  If  the  deliver 
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every  N  days  option  is  used,  the  next  delivery  is  schedule 
each  time  the  delivery  event  (or  subroutine)  is  executed. 

The  simulation  keeps  track  of  who  ordered  each  requisition 
so  when  it  arrives  it  can  be  counted  against  the  proper  com¬ 
ponent  under  repair.  It  also  keeps  track  of  the  progress  of 
the  components  so  that  when  repair  is  completed,  delay  costs 
can  be  assessed  to  the  proper  customer. 

This  simulation  model  accepts  an  arbitrary  number  of 
customers,  each  with  its  own  delay  cost  rate,  induction  per¬ 
iodicity,  and  demand  rate,  as  well  as  the  system  delivery 
cost  and  periodicity.  An  arbitrary  issue  delay  or  response 
time  function  can  also  be  specified. 

The  simulation  was  used  to  generate  points  on  the  deliv¬ 
ery  frequency-average  total  cost  curve,  with  random  number 
generator  seeds  being  reset  for  each  set  of  parameters  to 
reduce  variability  between  simulations.  The  simulation  was 
allowed  to  reach  steady  state  before  initializing  counters 
for  statistics  and  was  then  allowed  to  run  for  at  least  an 
additional  360  work  days.  The  simulation  was  based  on  24 
hour  work  days  and  ignored  the  effects  of  customers  not 
working  on  weekends  and  holidays. 
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A.  FIXED  DEMAND  HATE-RANDOM  ISSUE  DELAY  MODEL 

First  consider  the  case  where  there  is  a  random  issue 


processing  delay  for  each  requisition  submitted.  Since  much 
of  the  variation  in  delay  as  seen  by  the  customer  comes  from 
the  fact  that  all  required  material  may  not  be  available 
locally,  non-local  issues  must  be  considered.  This  enlarged 
system  is  illustrated  in  Figure  7. 
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Figure  7:  Enlarged  Customer-Stock  Point  Relationship 


To  construct  the  issue  delay  function  some  assumptions 
on  stock  point  effectiveness  and  system  responsiveness  were 
required.  It  was  assumed  that  the  local  stock  point  would 
fill  and  deliver  to  the  transportation  officer  55*  of  the 
requisition  submitted  in  2.5  days,  and  another  5*,  delayed 
for  some  unknown  reason,  would  be  filled  and  sent  to  local 
delivery  uniformly  throughout  the  next  4.5  days.  After 
delivery  to  the  transportation  officer  or  local  delivery, 
the  requisitions  would  be  delivered  to  the  customer  on  the 
next  scheduled  delivery.  The  2.5  day  local  issue  delay 
value  is  based  upon  some  requisition  processing  time  at  the 
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NARP  plus  standard  issue  processing  at  the  supply  center. 

The  Uniform  Material  Movement  and  Issue  Priority  System 
(UMNIPS)  standard  for  processing  issue  group  two  documents 
is  2  days.  This  can  be  improved  through  management  atten¬ 
tion  and  NSC  Oakland  has  made  it  a  policy  to  deliver  all 
maintenance  related  material  to  NARF  Alameda  within  one  day. 
The  60*  gross  effectiveness  at  the  local  stock  point  may 
seem  optimistic,  but  it  should  be  realized  that  data  col¬ 
lected  by  Hrabosky,  Owen,  and  Popp  [Ref.  6]  showed  that 
prior  to  consolidation,  NSC  Oakland  was  filling  36*  of  NARF 
Alameda  referrals.  This  36*  plus  whatever  Naval  Air  Station 
Alameda  was  filling  from  stocks  now  carried  by  NSC  Oakland 
may  give  the  60*  effectiveness  hypothesized. 

It  was  also  assumed  that  an  additional  25*  of  the 
requested  material  would  be  available  in  the  system  and 
would  be  shipped  by  non-local  means  directly  to  the  cus¬ 
tomer.  It  was  assumed  the  material  would  be  received  some¬ 
where  between  7  and  15  days  after  requisitioning.  It  was 
assumed  the  remaining  15*  of  the  items  required  would  be  out 
of  stock  and  the  backorder  and/or  procurement  process  would 
increase  delivery  time  to  the  customer  to  somewhere  uni¬ 
formly  distributed  between  15  and  45  days.  No  repair  parts 
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ware  assumed  to  have  a  leadtime  in  excess  of  45  days  (per¬ 
haps  another  optimistic  assumption)  . 

Because  non-local  deliveries  are  being  introduced  into 
the  system  and  because  only  the  local  delivery . system  is 
being  explicitly  modeled,  a  further  modification  to  the 
simulation  was  required.  For  parts  issued  non-locally, 
delivery  time  was  included  in  the  stated  leadtime  estimates 
while,  for  local  issues,  total  requisition  delay  is  the  sum 
of  issue  processing  time  and  the  time  to  make  the  delivery. 
As  a  consequence,  delay  costs  were  divided  into  two  compo¬ 
nents:  1)  those  caused  by  locally  issued  material  and  2) 
those  caused  by  non-local  issues. 

Because  requisitions  now  have  individual  lead  times,  the 
number  of  requisitions  submitted  has- now  become  a  factor  in 
delay  costs.  For  this  reason  each  production  line  supported 
has  a  new  parameter,  D,  which  is  the  number  of  requisitions 
submitted  per  component  repaired. 

&  two-customer  simulation  was  performed  with  parameter 
values  CT*100  dollars,  CD1*100  dollars,  11*7  days,  D1*14 
requisitions  per  component,  CD2*30  dollars,  Y2*3  days,  and 
02*6  requisitions  per  component.  Figure  8  shows  the  overall 
costs  and  delay  costs  contributed  by  local  and  non-local 
issues.  This  curve  is  being  considered  continuous  even 
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though  for  the  deterministic  case  it  was  not.  With 
probabilistic  variables  in  the  model,  the  perfect  phasing  of 
requisition  arrivals  are  no  longer  apparent  in  the  model. 
Moreover,  there  are  no  longer  the  wild  increases  and 
decrease  in  costs  noted  in  the  simulation  results.  This 
continuity  assumption  will  be  made  for  all  cost  curves  gen¬ 
erated  through  simulation  in  this  thesis. 

Several  interesting  facts  can  be  noted  from  Figure  8. 
First,  local  delivery  delay  costs  take  major  jumps  at  N 
values  of  11,  16,  and  20,  although  these  jumps  are  matched 
by  decreases  in  non-local  delay  costs.  These  jumps  are 
caused  by  the  discontinuities  in  the  issue  delay  probability 
distribution  function  and  are  believed  to  have  no  further 
significance. 

Hext,  the  delay  costs  have  driven  total  costs  much 
higher  than  in  the  previous  chapter.  Although  response  time 
for  local  issues  has  been  increased  2.5  days,  most  of  the 
delay  costs  are  now  coming  from  non-local  issues.  The  non¬ 
local  delay  costs  dominate  the  total  delay  costs  for  deliv¬ 
ery  schedules  of  10  days  or  less  resulting  in  a  much  flatter 
total  cost  curve  than  before.  It  is  only  when  locally 
delivered  material  begins  to  arrive  after  non-local  issues 
that  total  delivery  costs  begin  to  climb.  It  is  likely  that 
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COST  COST 


DELIVERY  PERIODICITY  (DRYS) 


QELIYERY  PERIODICITY  (DRYS) 


Figure  8:  Fixed  Deaand  R an doe  Issue  Delay  Case 
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the  local  stock  point  would  find  such  a  schedule  to  be 
"undesirable11  and  would  schedule  deliveries  sore  frequent 
than  is  optimal.  However,  because  the  cost  curve  is  so  flat 
the  stock  point  would  not  be  noticeably  affecting  the  total 
cost. 

It  should  also  be  noted  that  component  delay  is  now  a 
function  of  the  number  of  requisitions  submitted  per  compo¬ 
nent  because  of  the  assumed  gross  effectiveness  values. 
Pigure  9  is  a  graph  of  Customer  I's  (14  repair  parts  per 
component)  and  Customer  2's  (6  repair  parts  per  component) 
average  component  delays  in  days  versus  delivery  periodic¬ 
ity.  As  might  be  expected,  the  components  which  require 
fewest  parts  have  a  greater  sensitivity  to  delivery  sche¬ 
dules  because  they  are  more  likely  to  have  all  repair  parts 
available  locally. 

B.  RANDOM  COANTITT  DEMANDED  CASE 

Next  consider  the  case  where  the  number  of  requisitions 
submitted  per  component  is  random  and  the  issue  delay  func¬ 
tion  is  still  in  effect.  In  this  case  the  number  of  requi¬ 
sitions  per  component  is  described  by  a  probability  density 
function.  The  underlying  cause  of  a  requisition  is  a  failed 
repair  part  which  is  currently  installed  in  the  component 
under  repair.  Earlier  studies  [Ref.  3  and  4]  considered  the 
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Figure  9:  Customer  Delay  Comparison 


determination  of  the  need  for  a  repair  part  to  be  a  Ber¬ 
noulli  trial  where  the  repair  part  would  be  replaced  with 
probability  p.  For  multiple  like  repair  parts  the  sum  of 
Bernoulli  trials  with  a  common  p  forms  a  binomial  distribu¬ 
tion.  If  parts  are  not  alike,  then  the  p's  can  be  expected 
to  be  different  and  there  is  no  nice  distribution  for 
arbitrary  p. 
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For  comparison  purposes  the  same  two  custoaers  froa  the 
previous  section  were  modified  so  that  each  would  have  like 
coaponents  with  probability  of  repair  part  failure  of  0.5. 
By  fixing  the  failure  probability,  the  aean  demand  for  each 
customer  (D1  or  D2)  was  used  to  calculate  the  number  of  Ber¬ 
noulli  trials  or  installed  repair  parts  per  component. 
Figure  10  coapares  the  average  cost  curves  under  this 
modification  with  those  of  the  previous  section. 


OELIVERT  PERIODICITY  (CRTS) 


Figure  10:  Bandoa  Demand-Fixed  Demand  Cost  Comparison 
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As  can  be  seen,  costs  have  decreased  about  20  percent 
from  the  earlier  fixed  demand  case.  This  is  because  random 
demand  has  created  a  variation  in  the  number  of  requisitions 
submitted.  Since  the  binomial  distribution  is  symmetric  for 
p=0.5,  this  reduction  is  likely  due  to  the  increase  in  vari¬ 
ance  in  the  number  of  requisition  per  component.  Apparently 
the  benefit  from  one  fewer  requisition  exceeds  the  cost  of 
one  additional  requisition  per  component  repaired.  For  a 
distribution  with  decreasing  probabilities  in  the  tails  like 
the  binomial,  this  seems  logical.  As  the  number  of  requisi¬ 
tions  increases,  each  one  has  a  lower  probability  of  being 
the  critical  "last  item  received"  which  actually  determines 
delay  costs.  On  the  other  hand,  as  fewer  requisitions  are 
required,  the  probability  of  not  ordering  the  item  which 
would  have  determined  delay  (i.e.  the  probability  of  reduc¬ 
ing  delay  cost)  increases  at  an  increasing  rate,  when  N, 
the  delivery  periodicity,  becomes  large  this  argument  can 
break  down  since  being  in  the  lower  tail  of  the  distribution 
becomes  much  less  advantageous.  In  fact,  if  W  were  such 
that  a  local  issue  took  as  long  to  receive  as  a  non-local 
issue,  distribution  variance  should  make  no  difference.  It 
is  doubtful  a  stock  point  would  let  local  service  degrade  to 
this  level,  however. 
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A  probability  of  repair  part  failure  of  0.5  cannot 
always  be  assumed  and  the  p  value  does  have  significant 
effect  on  distribution  variance  and  shape.  To  investigate 
the  effects  of  a  varying  p  parameter,  a  siaulation  was  run 
with  3  customers  having  identical  parameters  of  CD*100,  T=7, 
and  D=6  and  only  the  binomial  distribution  p  values  were 
allowed  to  be  different.  Since  it  was  previously  shown  that 
components  which  required  more  repair  parts  had  higher  delay 
costs,  distribution  means  were  made  equal  to  6  by  varying  n, 
the  maximum  number  of  repair  parts  that  might  need  replac¬ 
ing,  along  with  the  paraaeter  p.  Customer  1  was  assigned  a 
p  value  of  0.1,  Customer  2  a  value  of  0.5,  and  Customer  3  a 
value  of  0.857.  Pigure  11  is  a  graph  of  the  average  compo¬ 
nent  delay  in  days  for  each  customer.  Included  in  the  graph 
are  the  delay  costs  experienced  in  the  p=J.O  or  determinis¬ 
tic  demand  case. 

The  customer  with  p  equal  0.5  generally  has  the  lowest 
delay,  the  one  with  p  of  0.  1  the  second  lowest,  and  the  one 
with  p  of  0.857  the  highest.  As  N  gets  large,  the  ranking 
is  not  so  clear,  however.  As  this  happens,  the  "long  lead 
time"  non- locally  issued  requisitions  actually  begin  arriv¬ 
ing  before  the  locally  issued  items.  It  appears  this  say  be 
becoming  a  problem  at  N*20. 
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Figure  11:  Delay  per  Component  For  Varying  Failure  Probability 


Por  N  less  than  20  there  are  at  least  two  forces  at 
work,  skewness  and  variance.  Pith  a  p  value  of  0.5,  the 
binoaial  distribution  is  syaaetric  about  its  mean  (has  zero 
skewness)  so  its  aean  is  also  the  median.  Thus  equal  num- 
bers  of  components  require  aore  than  the  aean  number  of 
requisitions  and  fewer  than  the  aean  nuaber.  As  discussed 
earlier  in  this  section,  a  few  sore  requisitions  increase 
delay  less  than  fewer  requisitions  reduce  it,  so  there  is  a 
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net  decrease  in  waiting  tiae.  As  p  decreases  from  0.5  sev¬ 
eral  things  happen.  First  n,  the  maximum  number  of  possible 
failed  parts  increases  because  the  mean  number  of  failures, 
np,  was  held  constant.  This  means  it  will  be  possible  for 
some  of  the  components  to  need  a  large  number  of  repair 
parts  (perhaps  n)  and  which,  as  a  consequence,  will  dominate 
delay  costs. 

Variance  is  also  affected  by  changing  n.  For  the  bino¬ 
mial  distribution  variance  is  np(l-p)  ,  and  since  np  is  being 
held  constant,  the  smaller  p,  the  larger  the  variance 
becomes.  Very  small  p  values  do  have  some  traits  which  tend 
to  increase  delay  costs.  The  distribution  does  become 
skewed  so  that  the  median  is  less  that  the  mean.  This  means 
that  the  decreases  from  the  mean  are  more  frequent  but  less 
in  magnitude.  Deviations  above  the  mean  are  infrequent  but 
are  quite  expensive.  These  are  the  inductions  which  domi¬ 
nate  costs  as  mentioned  above. 

Attempts  were  made  to  determine  where  the  decrease  costs 
from  increased  variance  were  overcome  by  the  effects  of 
higher  distributional  moments.  Simulations  were  run  with  p 
values  of  0.9,  0.67,  0.6,  0.4,  0.25,  and  0.2.  The  corres¬ 
ponding  n  values  were  7,  8,  9,  10,  15,  24,  and  30,  respec¬ 
tively.  The  differences  in  delay  costs  were  so  slight  that 
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no  strict  ordering  in  costs  could  be  obtained.  All  delay 
costs  sere  below  that  observed  in  the  deterministic  case, 
however. 

C.  THE  EFFECT  OF  VARIOUS  FAILURE  PROBABILITIES  IN  ONE 

COMPONENT 

In  the  previous  section  it  was  assuaed  that  all  parts 
within  a  component  had  the  same  probability  of  replacement 
(or  failure) ,  p.  It  was  also  hypothesized  that,  barring 
effects  from  higher  moments,  an  increase  variance  leads  to 
slightly  decreased  delay  costs.  By  examining  the  demand 
distribution  for  a  component  with  two  different  p  values, 
perhaps  a  statement  can  be  made  concerning  delay  cost 
estimates. 

Assume,  for  example,  a  component  had  2  classes  of 
required  repair  parts,  each  with  a  different  population 
(call  them  nl  and  n2)  and  a  different  probability  of  failure 
(pi  and  p2) .  Let  the  mean  of  the  approximate  distribution 
be  equal  to  the  sum  of  the  two  exact  binomial  distributions* 
Letting  n1+n2*n,  the  aggregate  demand  parameter  p  can  then 
be  defined  as 

p*  (nlpl  ♦  n2p2)  /n 

To  compare  variances,  the  sum  of  the  variances  of  the 
exact  distributions  should  be  compared  to  the  variance  of 
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the  approximate  distribution.  For  the  exact  distributions 
Var  D  (exact)  =  nlpl(l-pl)  ♦  n2p2(1-p2) 

=  nlpl  -  nlpl2  ♦  n2p2  -n2p22 
For  the  approximate  case  the  variance  is 
Tar  D  (approx)  =  np  (1-p) 

*  (ntpl  +  n2 p2)  (1-(n1p1  ♦  n2p2)/n) 

*  nlpl  ♦  n2p2  -  (nlpl  ♦  n2p2)  2/n 

Next,  set  the  difference  between  these  two  variances  equal 
to  a  constant  and  attempt  to  determine  the  sign  of  that 
constant. 

K  *  Tar  D  (approx)  -Tar  D  (exact) 

K  *  nlpl  ♦  n2p2  -  (nlpl  ♦  n2p2)  2/n  -  nlpl  -  n2p2 

♦  nlpl2  +  n2p22 

K  »  -(nlpl  ♦  n2p2)  2/n  ♦  nlpl*  ♦  n2p22 
nK  »  -  (nlpl  ♦  n2p2) 2  ♦  (nt  +n2)  (nlpl2  +  n2p22) 

»  -  n12p12  -  2n1n2p1p2  -  n22p22  ♦  n12p12 

*  n1n2p12  +  ntn2p22  ♦  n22p22 

*  n1n2  (pi*  -2p1p2  ♦  p22) 

*  n  1n2  (pi  -  p2)  2 

For  positive  nl  and  n2,  nK  and  thus  K  must  be  positive, 
indicating  the  variance  in  the  number  of  requisitions  sub¬ 
mitted  in  the  approximate  case  must  be  greater  than  in  the 
exact  case.  Although  the  above  argument  was  for  only  two 
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binomial  random  variables,  it  can  be  generalized  to  an 
arbitrary  number  of  values  for  p.  Thus  by  using  an 
appropriate  binomial  distribution,  variance  is  being 
understated  and,  ignoring  the  effects  of  higher  moments, 
delay  costs  are  being  overestimated. 

D.  EXAMINATION  OF  THE  "SHIP  EVERY  K  REQUISITIONS” 

PHILOSOPHY 

As  stated  earlier,  a  study  by  Davidson  [Ref.  3]  shoved 
little  difference  in  the  optimal  costs  for  the  local  deliv¬ 
ery  options  listed  at  the  beginning  of  Chapter  3.  To  verify 
this  in  the  multi-customer,  multi-item  inventory  case,  simu¬ 
lations  vere  run  to  compare  the  "Ship  every  N  days"  strategy 
to  the  "Ship  every  K  requisitions"  philosophy. 

In  making  comparisons  between  these  plans,  some  sort  of 
equivalency  must  be  developed.  Comparing  a  plan  where  K=10 
with  an  N=2  may  give  one  result  when  the  system  is  deliver¬ 
ing  roughly  5  repair  parts  per  day  and  quite  another  if  on 
the  average  50  repair  parts  per  day  are  being  shipped.  For 
this  reason  it  was  decided  to  compare  plans  where  the  mean 
numbers  of  parts  per  delivery  were  approximately  equal. 

Under  the  deliver-every-K-requisitions  option,  obviously  the 
load  is  always  K  requisitions.  For  the  deliver-every-N-days 
case,  the  mean  delivery  load  is  the  average  daily  demand 
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tiaes  the  proportion  of  requisitions  shipped  via  local 
delivery  (0.6  with  the  previously  defined  issue  delay  func¬ 
tion)  times  the  number  of  days  between  deliveries.  Express¬ 
ing  this  aatheaat ically , 

Average  Delivery  *  Average  Daily  •  (0.6)  ♦  (N)  .  (4.1) 

Load  Deiand 

Por  each  custoaer  the  average  daily  demand  would  be  the 

average  number  of  requisitions  per  repair  divided  by  the 

period  between  repairs.  Summing  this  for  average  daily 

demand  for  the  two-customers  case. 

Average  Daily  *  n  1p1  ♦  n2p2  .  (4.2) 

Demand  Y  i  1 2 

Combining  equation  (4.1)  and  equation  (4.2), 

Average  Delivery  =  (0.6)  •(**)•  fnl p  1  ♦  n2p2\  .  (4.3) 

Load  ITi  t'2  I 

Costs  for  delivery  plans  with  equivalent  average  load  values 

can  now  be  compared. 

Two-customers  simulations  were  run  with  parameters 
p1*0. 1  ,  p2  =  0. 1 ,  n 1=175,  n2=75,  Y1  =  7,  12=3  ,  CT=100,  CD1=100, 
and  CD2=30.  Dsing  equation  (4.3)  it  can  be  seen  the  average 
delivery  load  should  be  3ti,  or  a  plan  with  H=3  should  be 
compared  with  a  plan  where  K=9. 

In  Pigure  12  average  total  costs  were  plotted  against 
the  average  nuaber  of  components  per  delivery  for  the  two 
plans.  These  cost  curves  are  nearly  coincident  and  thus  it 
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Figure  12:  Cost  Comparison  for  Two  Delivery  Strategies 

appears  the  two  delivery  strategies  are  equivalent.  It 
should  be  reaeabered,  however,  that  it  was  assuaed  local 
delivery  was  not  capacity  constrained.  If,  in  the  delivery- 
every-N-days  aoda,  aaterial  was  not  delivered  because  of  a 
capacity  constraint,  then  delay  costs  would  be  higher.  For 
the  deliver-every-K-requisition  strategy  such  a  problea 
could  not  exist  since  delivery  capacity  aust  be  at  least  K 
requisitions  for  feasibility. 


V.  amis  £M2  5I2C5M2  hi  IBS  ifiS-ssimi  SIIS 


The  last  problea  to  be  considered  is  whether  a  separate 
warehouse  facility  should  be  maintained  at  the  industrial 
repair  facility  to  support  operations.  Although  the  model 
does  not  determine  what  should  be  stocked  and  where,  the 
costs  of  various  alternatives  can  be  analyzed.  For  example, 
if  establishing  a  local  warehouse  will  decrease  delay  costs 
for  one  NARP  production  line  much  more  than  another,  perhaps 
the  local  warehouse  should  concentrate  on  carrying  stock  for 
the  line  which  derives  the  greater  benefit.  First,  though, 
a  more  basic  question  must  be  asked. 

A.  FACTORS  APFECTI NG  THE  REMOTE  WAREHOUSE  DECISION 

Hany  factors  contribute  to  the  decision  of  whether  or 
not  to  establish  a  customer-sited  warehouse.  A  review  of 
some  of  these  factors  will  place  the  delay  cost  problea  in 
perspective. 

First,  the  overhead  of  maintaining  a  separate,  remote  • 
sited  warehouse  must  be  considered.  Extra  material  and  per¬ 
sonnel  are  likely  to  be  required.  If  automation  in  handling 
materials  at  the  main  warehouse  has  made  it  more  efficient. 
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the  increased  costs  of  daily  processing  at  the  remote  ware¬ 
house  eust  be  included. 

Next,  the  source  of  aaterial  scored  in  the  remote  ware¬ 
house  aust  be  considered.  If  material  is  received  primarily 
from  off-base  and  is  routed  through  a  centralized  receiving 
at  the  main  supply  center,  handling  the  material  at  both  the 
main  center  and  the  remote  warehouse  can  incur  significant 
extra  cost.  If,  on  the  other  hani,  aaterial  represents  com¬ 
ponents  which  have  been  made  ready-for- issue  by  the  indus¬ 
trial  facility  being  supported  and  which  are  being  returned 
to  the  system,  significant  savings  can  result  by  stocking 
the  material  at  the  remote  site.  Ihis  is  especially  true  if 
the  aaterial  is  issued  to  another  customer  at  the  remote 
site,  such  as  another  N&RP  production  line. 

The  speed  of  stock  record  take-up  is  another  important 
factor,  although  costs  are  diff  icult  to  quantify.  By  avoid 
transshipment  of  repaired  material  to  the  main  supply  cen¬ 
ter,  stock  records  can  be  updated  sooner  and,  if  require¬ 
ments  for  the  repaired  component  exist,  the  issue  can  be 
made  more  rapidly. 

But  perhaps  the  most  obvious  benefit  of  stockage  of 
material  at  the  customers*  site  is  t^he  decrease  in 
requsition  waiting  time.  By  modifying  the  issue  delay  time 
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function  introduced  in  the  last  chapter,  delay  cost  inpacts 
can  be  estimated. 

8.  SIMULATING  THB  SYSTEM  WITH  A  LOCAL  HABBHOOSE 

In  simulating  a  system  vith  a  local  warehouse,  only  the 
issue  delay  function  needs  to  be  aodified.  It  was  assumed 
that  if  such  a  warehouse  existed  it  would  fill  of  the 
demands  submitted  by  the  co-located  industrial  customers. 
This  gross  effectiveness  is  ju3t  a  rough,  perhaps  pessimis¬ 
tic  guess  at  what  might  be  obtained  by  a  standard,  demand 
based,  stocking  policy.  By  making  issues  locally,  material 
would  not  have  to  enter  the  supply  center's  local  delivery 
system  and,  it  was  assumed,  would  be  available  to  the  cus¬ 
tomer  in  exactly  one  day.  The  issue  delay  function  was 
modified  accordingly  and  the  simulation  was  run  for  four 
co-located  customers  using  the  policy  "ship  every  N  days". 
Four  customers  were  chosen  to  provide  a  spread  in  customer 
parameter  values.  The  number  of  demands  per  component  for 
each  customer  were  binomially  distributed  with  p=0.1  and  all 
customers  were  assigned  a  delay  cost  rate  (CD)  of  100.  Y 
(the  time  between  inductions)  and  D  (number  of  requisitions 
per  induction)  were  equal  for  each  customer  but  were  diffe¬ 
rent  for  each  of  the  four,  being  18,  12,  6,  and  3,  respec¬ 
tively.  These  values  of  Y  and  D  allowed  each  customer  to 
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have  an  average  demand  rate  of  one  requisition  per  day,  yet 
provide  a  spread  in  the  average  number  of  requisitions  per 
coaponent.  Figure  13  is  a  comparison  of  the  local  warehouse 


DELIVERY  PERIODICITY  (DAYS) 

Figure  13:  Local  vs  Non-Local  Stocking 


siaulation  versus  the  no  local  warehouse  simulation.  As  can 
be  seen,  costs  were  roughly  1  to  '4  percent  lower  for  the 
local  warehouse  case  and  only  1  to  2  percent  lower  for 
values  of  N  less  than  ten. 


Next  it  could  be  asked  which  of  the  above  four  customers 


benefitted  the  nost  from  the  local  warehouse.  If  Custoaer  4 
(3  demands  per  repair  on  the  average)  showed  a  significant 
delay  cost  reduction,  perhaps  more  of  his  material  should  be 
stocked  in  the  local  warehouse  even  at  a  cost  of  having  less 
material  for  the  other  customers.  Maybe  material  should  be 
stocked  so  that  all  of  his  supply  canter  issues  should  be 
made  from  the  local  warehouse  while  only  a  few  issues  are 
made  locally  for  the  other  customers.  Of  course,  stocking 
to  a  higher  effectiveness  usually  requires  higher  and  higher 
investment  per  incremental  issue,  and  perhaps  a  cost-benefit 
analysis  is  appropriate. 

Using  the  same  simulation  model  as  for  the  previous  fig¬ 
ure,  individual  custoaer  average  component  delays  were  cal¬ 
culated  for  each  warehousing  plan.  For  both  Customer  1  and 
Custoaer  2,  average  component  delay  was  the  same  with  or 
without  the  remote  warehouse  for  all  values  of  N  between  2 
and  20  days.  Figure  14  shows  the  graphs  of  average  compo¬ 
nent  delay  for  the  other  two  customers. 

It  appears  Custoaer  4,  the  custoaer  who  on  the  average 
only  required  3  repair  parts  per  component  repaired,  would 
benefit  most  from  a  co-located  warehouse.  Customer  3  (6 
demands  per  repair)  would  also  benefit  some,  though  it 
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figure  14:  Customer  Delay  Cost  Coaparisons 
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appears  as  the  average  number  of  requisitions  per  component 
increase,  the  benefit  derived  from  a  local  warehouse 
decreases.  As  might  be  expected,  as  deliveries  become  more 
freguent,  local  warehouse  benefits  also  decrease.  Thus 
given  similar  delay  cost  rates  and  demand  rates,  if  a  local 
warehouse  has  been  established,  delay  costs  can  be  reduced 
by  targeting  stocked  material  to  the  customer  who  require 
the  fewest  repair  parts  per  component  repaired. 

C.  THE  EPFECTS  OF  IMPROVED  EFFECTIVENESS 

Throughout  earlier  analyses  it  has  been  assumed  the  sup¬ 
ply  center  has  been  limited  to  60*  point  of  entry  (POE) 
effectiveness.  Shat  would  happen  if,  by  studying  past  fai¬ 
lure  data  and  possibly  making  increases  in  range  and  depth, 
effectiveness  could  be  increased?  Assume,  for  example,  the 
supply  center  could  fill  75*  of  the  NARF  requisitions  in  2.5 
days  and  an  additional  5*  in  the  next  five  days.  If  the 
remaining  20*  of  the  requisitions  were  split  evenly  between 
system  issues  (7  to  15  days  from  requisition  date  until 
receipt  by  customer)  and  backorders  (15  to  45  days  until 
receipt) ,  a  new  issue  delay  function  is  defined.  Using  the 
same  four  customers  as  in  the  simulations  used  for  Figure 
13,  cost  curves  were  generated  for  this  new  issue  delay 
function.  The  new  curves  are  shown  in  Figure  15.  The 
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Figure  15:  Local  vs  Non-Local  Stocking-Enhanced  Effectiveness 

highest  curve  is  for  the  case  of  no  remote,  customer-sited 
warehouse.  The  middle  curve  is  for  a  remote  warehouse  which 
makes  two  thirds  of  the  supply  center’s  issue  to  the  NARF, 
the  same  percentage  as  was  considered  earlier.  The  lowest 
curve  represents  a  remote  warehouse  which  is  able  to  make 
75*  of  the  supply  center's  issues  to  the  NA5F. 

As  can  be  seen  by  comparing  Figures  13  and  15,  with  the 
higher  effectiveness  total  daily  costs  for  N  between  2  and  6 
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have  dropped  nearly  22%  for  the  no  local  warehouse  case. 
When  it  is  assuaed  the  local  warehouse  Bakes  two  thirds  of 
the  HSC’s  issues  as  before,  the  percent  cost  reduction  is 
even  slightly  greater.  Further  iaprovement  is  possible, 
however,  by  assuming  the  local  warehouse  can  sake  three 
quarters  of  the  NSC's  issues  to  the  NARF.  This  is  shown  by 
the  lowest  curve  in  Figure  15.  Since  more  issues  are  made 
locally,  it  is  more  likely  all  material  is  available  and 
rapid  local  response  can  be  converted  into  lower  delay 
costs.  Thus  it  appears  that  success  feeds  upon  itself  and 
those  activities  with  the  highest  effectiveness  can  benefit 
the  most  from  a  remote  warehouse. 
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71.  SUMMARY  AND  CONCLUSIONS 


Unfortunately  this  paper  was  unable  to  find  a  simple 
solution  or  algorithm  for  the  optimal  delivery  or  siting  of 
repair  parts.  When  the  local  delivery  problem  is  general¬ 
ized  to  a  multi-customer,  multi-repair  part  inventory 
environment,  the  number  of  different  parameters  becomes 
significant  and  even  in  the  relative  simple  deterministic 
case  the  cost  function  lacks  continuity  as  well  as 
convexity.  Nevertheless,  there  is  an  upper  bound  function 
which  can  be  optimized.  This  bounding  function  has  a  well 
defined  minimum  which  is  similar  in  form  to  the  cost 
function  in  the  Hadley  and  Whitin  Deterministic  Lot 
Inventory  Model.  It  was  also  noted  that  the  number  of 
requisitions  submitted  per  component  repaired  had  no  effect 
on  costs. 

In  generalizing  to  the  multi-item  inventory,  which 
allowed  more  than  one  requisition  per  component  repaired,  a 
key  assumption  was  made  concerning  the  assessment  of  delay 
costs.  By  allowing  delay  costs  to  accumulate  at  a  constant 
rate  until  all  ordered  parts  were  received,  much  more  empha¬ 
sis  was  placed  on  requisitions  with  the  slowest  delivery 
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tines.  When  a  function  which  simulates  system  response 
times  as  veil  as  local  issue  processing  times  was 
constructed#  these  slow  requisitions  were  generally  issues 
outside  the  local  system  and  thus  not  a  function  of  the 
local  delivery  schedule.  As  a  consequence#  optimal  costs 
became  very  insensitive  to  delivery  schedules.  At  the  other 
end  of  the  spectrum,  delay  costs  became  sensitive  to  the 
number  of  requisitions  ordered  per  component#  particularly 
when  that  number  was  small.  Those  components  requiring  few 
parts  could  more  often  have  all  requirements  filled  at  the 
requisition  point  of  entry  (POE)  and  thus  experience  minimal 
delay.  An  increased  POE  effectiveness  would  also  provide  a 
similar  decrease  in  delay  cost. 

Delay  costs  were  also  decreased  when  variability  was 
allowed  in  the  number  of  requisitions  per  component.  Vari¬ 
ance,  though#  was  not  the  only  distributional  moment  which 
affected  delay  costs,  for  costs  also  seemed  to  increase  as 
the  p  decreased  below  0.5  using  the  binomial  demand  distri¬ 
bution.  More  study  should  be  conducted  in  this  area. 

Lastly,  the  warehousing  of  material  at  remote  sites  was 
considered.  As  modeled#  a  warehouse  located  at  the  custom¬ 
ers*  site  had  little  impact  unless  exceedingly  few  parts 
were  required  per  component.  The  model  assumed  only  40%  of 
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all  demands  could  be  filled  from  the  local  warehouse  and  a 


greater  effectiveness  could  give  greater  delay  cost 
reductions.  It  appears  high  effectiveness  and  targeting 
material  towards  customers  who  require  only  a  few  repair 
parts  per  component  is  essential  to  deriving  full  benefits 
from  remotely  located  warehouses. 

With  respect  to  customer  response  time  the  following 
conclusions  can  be  made.  First,  this  study  shows  that  non¬ 
local  deliveries  and  POE  effectiveness  are  usually  the  lim¬ 
iting  factors  in  delay  costs.  Although  most  issues  for 
local  customers  will  not  reduce  waiting  time,  many  issues  to 
non-local  customers  may  be  critical  "last  part  required"  and 
thus  reduce  system  delay  costs.  This  means  that  when  an 
activity  such  as  NSC  Oakland  invests  in  equipment  which 
reduces  response  time,  not  only  are  delay  costs  reduced  at 
local  customers  such  as  NARF  Alameda,  but  there  also  may  be 
reductions  at  other  major  customers  such  as  Ship  Repair 
Facility,  Subic  Bay,  Philippines,  or  Ship  Repair  Facility, 
Yokosuka,  Japan.  Expeditious  deliveries  to  fleet  units 
located  at  the  industrial  site  are  important,  since  the  lack 
of  repair  parts  may  be  directly  affecting  fleet  readiness. 

Next,  as  might  be  expected,  the  more  requisitions  ord¬ 
ered,  the  greater  the  delay  cost.  Thus  if  many  required 
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iteos  are  stocked  as  planned  re guirements,  pre-expended  bin, 
or  in  repair  kits,  the  total  number  of  requisitions 
submitted  to  the  POE  at  the  time  of  repair  can  be  decreased, 
decreasing  delay  costs. 

By  their  nature,  delay  costs  are  somewhat  nebulous  and 
the  results  of  this  study  could  be  challenged  on  those 
grounds.  Perhaps  having  some  repair  parts  rapidly  available 
would  decrease  delay  costs.  Perhaps  an  upper  bound  (or  time 
standards)  on  supply  response  time  is  appropriate  such  that 
delay  costs  would  only  be  assessed  when  this  time  is 
exceeded.  Unfortunat ely,  time  standards  are  now  dictated  by 
the  system  rather  than  by  individual  repair  processes.  The 
Jfaval  Avaiation  Maintenance  Plan  (N&MP)  says  only  that  issue 
group  one  material  must  be  delivered  within  an  hour  and 
issue  group  two  and  three  in  two  hours  and  twenty  four  hours 
respectively,  regardless  of  the  repair  process.  Lastly, 
perhaps  delay  costs  are  not  time  dependent  and  only  a  fixed 
charge  should  be  assessed  if  time  standards  are  not  met. 

More  investigation  on  the  nature  of  delay  costs  appears  in 
order. 

This  study  also  assumed  all  requisitions  were  treated 
equally  by  the  system.  Ihere  were  no  issue  priorities,  pre¬ 
mium  transportation,  or  material  expediters.  Expediting 
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critical  delay-causing  requisitions  would  be  a  particularly 
effective  way  of  decreasing  delay  costs  in  this  model.  This 
could  ideally  be  done  through  computerized  requisition  sub¬ 
mission,  follow-up  and  monitoring  programs.  Only  through 
good  local  requisition  processing  and  expedited  system  sup¬ 
port  can  industrial  facilities  keep  the  depot  turnaround 
time  to  a  minimum  and  operational  availability  at  a  maximum. 
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Consider  the  component  delay  problem  as  illustrated  in 
Figure  2  of  Chapter  3.  As  long  as  both  H  (the  number  of 
days  between  deliveries)  and  Y  (the  number  of  days  between 
inductions  or  equivalently  between  requisition  submissions) 
are  constant  and  H/Y  is  a  rational  number,  component  delays 
are  cyclic.  Moreover,  the  average  delay  per  component  over 
the  cycle  can  be  calculated. 

Theorem  1:  If  N/Y  is  a  rational  number  and  N  and  Y  are 
constants,  the  values  for  delay  cost  will  be  cyclic  over 
time. 

Proof:  The  figure  below  shows  a  timeline  of  two  compo¬ 
nent  repairs  where  D(0)  and  D(1)  are  the  delay  tines  for  two 
consecutive  components. 

t— D  (0)— H  T  1*-D  (1 )  — H 

0  N  (j  +  1)H 

Ojder  Order 

Figure  16:  Repair  Timeline 

The  key  points  to  notice  in  the  figure  is  that  deliveries 
are  an  integer  times  N  days  apart,  inductions  are  Y  days 
apart,  and  delays  are  measured  from  an  induction  to  the  next 
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delivery.  Proa  the  figure 

C  ( 1)  ♦  Y  =  D  (0)  ♦  jH  for  some  integer  j ,  or 
D(1)  =D(0)  ♦  j  N  -  Y  . 

For  the  i  th  component  inducted,  this  expression  becomes 
E(i)  =  D  (0)  ♦  jN  -iY  , 

where  j  is  an  integer  and  is  chosen  such  that  D(i)  is  the 
smallest  positive  number  possible.  But  since  N/Y  is 
rational, 

=  41_  for  some  integers  u  and  v. 

If  this  is  the  case,  uY  *  vN,  or  D(k)  =  D  (0)  since  j  can  be 
chosen  to  be  v.  Thus  delay  values  are  cyclic  every  u  induc¬ 
tions  and  the  cycle  length  is  uY  days. 

Although  as  drawn  it  appears  N  is  being  restricted  to  a 
value  less  than  Y,  this  is  not  necessarily  the  case.  If 
J=*1,  as  long  as  D  (0)  is  less  than  Y  (  as  it  must  be  for  some 
component  whose  last  repair  part  is  delivered  at  time  N) , 
the  above  and  below  arguments  hold,  although  the  figure  may 
not  be  to  scale.  D(1)  nay  accept  values  greater  than  Y. 

Theorem  2:  If  delay  costs  are  cyclic,  the  average  delay 
over  the  cycle  in  days  is 

D  *  N  -  L  ♦  g 

where  L  is  the  largest  real  number  such  that  N/L  and  Y/L  are 
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both  integer  and  g  is  soae  constant  between  0  and  L. 

Proof:  Since  delay  costs  are  cyclic,  let  D(0>  be  an 
arbitrary  delay  observed  in  a  cycle.  Subtracting  an  cons¬ 
tant  g  from  each  observed  delay  the  recursion  relation  from 
Theoren  1  above  becomes 

D (1)  -  g  =  D(0)  -  g  -  r  *  jN  or 

0(1)  -  g  =  n  ■  |d  ( o )  -  a  -  r  *  m 

Since  deliveries  occur  every  N  days  and  orders  are  filled 
iniediately,  the  aaxiaua  delay  will  be  no  aore  than  N  days. 
If  that  is  the  case,  the  expression  in  parentheses  aust  only 
assuae  values  less  than  1.  If  the  calculated  delay,  D(1)-g, 
is  greater  than  V ,  the  required  part  would  have  been  deliv¬ 
ered  with  an  earlier  delivery  N  or  soae  multiple  of  N  days 
earlier.  This  leaves  only  the  fractional  part  of  the  above 
expession  in  parentheses  as  delay.  The  expression  can  then 
be  rewritten 

D  (1 )  -  q  *  H  Fractional  Part  (D  (0 )  -  a  -  Y  ♦  dm 

N 

Mext  divide  both  the  numerator  and  the  denominator  of  the 
fraction  and  both  sides  of  the  equation  by  L,  where  L  is  the 
largest  real  nuaber  such  that  and  l/L  are  both  integer. 
The  expression  then  becomes 
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( -  g  =  N  Fractional  Part 


SHU  ?  a  -  i  ♦  in 


N/L 


*  (PiO.LL-  q 


mod  N 


(A.  1) 


‘  S'-f) 

where  sod  is  the  modulus  function.  Using  this  latest 
expression#  first  consider  the  term 


Since  j,  N/L,  and  T/L  are  all  integer,  this  expression,  call 
it  c,  must  be  an  integer.  Moreover,  since  this  tern  appears 
within  the  parentheses,  j  can  be  adjusted  so  that  c  accepts 
values  between  0  and  N/L  without  affecting  the  equation. 

Next  define  the  general  expression 
D  (i|^  -q  as  X  (i)  , 

which  can  be  used  to  describe  both  the  left  hand  side  of 
equation  (A.1)  and  the  first  tern  in  the  mod  expression. 
First  consider  X(0),  the  tera  on  the  left  side  of  the  sod 
expression.  By  chocsing  soae  q  between  0  and  L,  X (0)  can  be 
aade  an  integer.  Moreover,  if  this  g  were  determined  when 
D (0)  was  the  smallest  delay,  all  other  X(i)  will  be  positive 
even  with  this  q  subtracted  from  D (i) .  In  any  case,  by 
defining  N/L  as  a  positive  integer  a,  the  mod  expression 
becoaes 

X  ( 1)  ■  (X  (0)  *  c)  mod  a 
or  aore  generally. 
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X  (i+1)  =  (X  (i)  *  c)  iod  i  . 

This  is  a  special  condition  of  a  linear  congruential  nuaber 
generator,  which,  as  discussed  by  Knuth  [Ref.  7],  is  of  the 
fora 

X(n+1)  *  (aX(n)  ♦  c)  mod  a  , 

where  X(0)  ,  a,  and  c  are  non-negative  integers  and  a  is  an 
integer  greater  than  X(0),  a,  or  c.  These  generators  are 
said  to  be  full  cycle,  or  accept  integer  values  from  0  to 
a-1,  if  the  following  conditions  are  met: 

1.  c  has  no  priae  factors  in  coaaon  with  a. 

2.  a  aod  y  =  1  for  all  y  which  are  priae  factors  of  a. 

3.  a  aod  4  *  1  if  4  is  a  factor  of  a. 

For  the  component  delay  case  Conditions  2  and  3  are  aet 
easily  since  the  paraaetsr  a  has  value  1,  and  thus  has  no 
integral  factors  other  than  1.  To  check  Condition  1  it  aust 
be  shown  that  N/L  and ( { jN/L) - (Y/L))  have  no  coaaon  priae 
factors.  First  assuae  such  a  factor  exists  (  call  it  z) . 

For  z  to  be  a  factor  of  N/L,  N/zL  aust  be  an  integer.  Since 

N/zL  and  j  are  integers,  jN/zL  aust  be  an  integer.  If  z  is 

a  factor  of  ( (jN/L)  -  (Y/L) ) ,  then  ( ( jN/L)  -  (T/L) ) /z  aust  be  an 
integer  or  (jn/zL)- (Y/zL)  aust  be  integer.  But  it  has 
already  been  shown  that  jN/zL  is  integer  so  Y/zL  aust  be 
integer.  But  if  Y/zL  and  N/zL  are  both  Integer  for  a  z 
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greater  than  1,  then  L  was  not  chosen  properly  (it  was  not 
the  largest  real  number  such  that  N/L  and  Y/L  are  both 
integers).  Thus  ( (jH/L) - (Y/L) )  can  have  no  common  factors 
with  H/L  and  Condition  1  aust  hold.  The  variables  X(i)  must 
assuae  values  0,1 ,...,m-1  or  have  an  average  value  of 
(a-1)/2.  Converting  the  X  variables  bach  to  0  and  using  the 
fact  that  a  is  N/L, 


X  *  D  9  =  O-^-l 


D  -  q  =  (ML*  L) 


D 


♦  9 


which  is  the  average  delay  cost  per  component  over  the 
cycle.  Since  delay  is  being  minimized  in  this  thesis,  it  is 
assumed  initial  conditions  will  be  established  such  that  g 


is  egual  to  zero. 
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PREAMBLE 

"OSES  ISSDELT  FUNCTION 

EVENT  NOTICES  INCLUDE  DELIVERY,  START. STAT,  END. SIM 
EVERI  RECEIPT  HAS  A  GDOC 
EVERY  INDUCTION  HAS  A  SNCUST 
EVERY  ISSUE  HAS  AN  GITEM 

PRIORITY  ORDER  IS  START. STAT,  INDUCTION,  ISSUE,  DELIVERY, 
RECEIPT,  AND  END. SIM 

DEFINE  GDOC,  GNCUST  AND  GITEM  AS  INTEGER  VARIABLES 
"  G  PREFIX  ON  VARIABLES  TO  DENOTE  GLOBAL  VARIABLE 
TEMPORARY  ENTITIES 

EVERY  COMPONENT  HAS  A  NPARTS,  A  S DELAY  AND  A  LINE,  AND 
BELONGS  TO  THE  REPAIR 

"  REPAIR  IS  THE  SET  OF  ALL  ITEMS  UNDER  REPAIR. 

DEFINE  NPARTS  AS  AN  INTEGER  VARIABLE 
EVERY  REQN  HAS  A  NHA  AND  MAY  BELONG  TO  THE  QUEUE 
"  QUEUE  IS  SET  OF  ALL  REQUISITIONS  A  WAITING  DELIVERY 
DEFINE  LINE  AND  NHA  AS  INTEGER  VARIABLES 
PERMANENT  ENTITIES 

EVERY  CUSTOMER  HAS  AN  INDMIN,  AN  INDMAX,  A  PMEAN, 

A  NBRRFI,  A  COSTRATE.  A  DLY1COST,  AND  A  DLY2COST 
"  DLY1COST  FOB  LOCALLY  DELIVERED  ITEMS,  A  DLY2COST  FOR  OTHERS 
THE  SYSTEM  OWNS  THE  QUEUE  AND  THE  REPAIR,  AND  HAS  AN  ISSDELAV 
RANDOM  LINEAR  VARIABLE 

DEFINE  IS3DELAY  AS  A  REAL,  STREAM  3  VARIABLE 
DEFINE  I.  PMEAN,  K,  NBRRFI,  TLOAD,  FLAG  AND  KK  AS  INTEGER 
VARIABLES 

DEFINE  TRANSCOST,  SHIPCOST,  STARTSTAT,  ENDSIM,  AND  N  AS 
VARIABLES 

TALLY  SDLY1COST  AS  THE  SUM  OF  DLY1COST 
TALLY  SDLY2COST  AS  THE  SUM  OF  DLY2COST 

TALLY  MN.N.QU  AS  THE  MEAN  AND  VAR.N.QU  AS  THE  VARIANCE 
OF  TLOAD 

END 


MAIN 

DEFINE  SS  AS  A  1- DIMENSIONAL  ARRAY 
RESERVE  SS{>  AS  4 

DEFINE  J,  NIT  AS  INTEGER  VARIABLES 
READ  NIT 

READ  N. CUSTOMER,  SHIPCOST 

PRINT  2  LINE  WITH  N. CUSTOMER,  SHIPCOST  AS  FOLLOWS 
N. CUSTOMERS  SHIPPING  COSTS 

***  ****.** 

READ  ISSDELAY 

CREATE  EVERY  CUSTOMER 

PRINT  1  LINE  AS  FOLLOWS 

CUSTOMER  INDMIN  INDMAX  PMEAN  COSTRAT] 

FOR  EACH  CUSTOMER,  DO 

READ  INDMIN  (CUSTOMER)  ,  IN DMAX  (CUSTOMER)  ,  PMEAN  (( 


COSTRATE 


INDMAX (CUSTOMER)  ,  PMEAN (CUSTOMER)  , 


COSTRATE (CUSTOMER) 

PRINT  1  LINE  WITH  CUSTOMER. INDMIN (CUSTOMER) , INDMAX (CUSTOMER) 
PMEAN  (CUSTOMER)  ,  COS  TR ATE  (CUSTOMER)  AS  FOLLOWS 

****  <***.*#  'A***.**  ****.**  ****.** 


LOOP 

FOR  I«1  TO  4,  DO 
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I  • 


SAVE  ALL  RANDOM  NUMBER  SEEDS 
LET  SS  (I)  “SEED.  V  ( I) 

LOOP 

POR  J=1  TO  NIT,  DO 

••  RUN  AS  MANY  ITERATIONS  AS  DESIRED 
READ  N.  K.  STARTSTAT,  ENDSIM 
POR  1=1  TO  4.  DO 
••  RESET  SEEDS  EACH  ITERATION 
LET  SEED.  V(I)  =SS(I) 

LOOP 

SKIP  2  LINES 
POR  EACH  CUSTOMER,  DO 
CREATE  A  COMPONENT 
PILE  THE  COMPONENT  IN  REPAIR 
LET  LINE  (COMPONENT)  =CUSTOMER 
LET  Q—  *1 

"  USING*BINOMIAL  WITH  P=.1.  PMEAN  IS  NUMBER  OF  TRIALS. 

LET  NPARTS  (COMPONENT)  =  B IN OMI A L.  F  (PMEAN  (CUSTOMER)  ,Q,4) 

LET  INDUCT=UNIFORM. F  (INDMINJCUSTOME R),  INDMAX (CUSTOMER) , 2) 
"  THIS  ALLOWS  VARIABILITY  POR  TIME  BETWEEN  INDUCTIONS. 

"  IN  THESIS  MAX  VALUE  AND  MIN  VALUE  WAS  ALWAYS  THE  SAME. 
LET  SDELAY(COMPONENT[=0 
POR  1=1  TO  NPARTS  (COMPONENT)  ,  DO 
CREATE  A  REQN 
LET  NHA  (REQN)  =COMPONENT 
LET  IWAIT=ISSDELAY 
LET  WAIT  =IWAIT 

IF  IWAIT  LE  7  AND  IWAIT  GE  2 . 5 
IF  IWAIT  IS  GT  7,  THE  ISSUE  IS  NON-LOCAL.  IF  IT  IS 
LT  2.5.  ISSUE  IS  FROM  LOCAL  WAREHOUSE.  IN  EITHER  CASE 
LOCAL  DELIVERY  SYSTEM  IS  NOT  USED. 

SCHEDULE  AN  ISSUE  GIVEN  REQN  IN  WAIT  DAYS 
ELSE  SCHEDULE  A  RECEIPT  GIVEN  REQN  IN  WAIT  DAYS 
REGARDLESS 
LOOP 

SCHEDULE  AN  INDUCTION  GIVEN  CUSTOMER  IN  INDUCT/2  DAYS 
LOOP 

SCHEDULE  AN  END. SIM  AT  ENDSIM 
SCHEDULE  A  START. STAT  AT  STARTSTAT 
LET  KK=  K 
LET  K=9999 

* »  THIS  K  VALUE  SERVES  AS  FLAG  THAT  DELIVERY  EVERY  N 
»•  OPTION  IS  IN  USE.  KK  SAVES  K  VALUE. 

SCHEDULE  A  DELIVERY  IN  N  DAYS 
START  SIMULATION 

"  REPEAT  FOR  DELIVER  EVERY  K  REQN  OPTION. 

LET  TIME. V=0 
IF  KK  NE  0 

FOR  1=1  TO  4,  DO 
LET  SEED.V  ( I)  *SS  (I) 

LOOP 

FOR  EACH  CUSTOMER,  DO 
CREATE  A  COMPONENT 
FILE  THE  COMPONENT  IN  REPAIR 
LET  LINE  (COMPONENT)  “CUSTOMER 
LET  Q=. 1 


i  • 


•  t 
t  • 


DAY 


LET  NPARTS  (COMPONENT)  “BINOMIAL.  F  (PMEAN  (CUSTOMER  1,0,4) 

LET  INDUCT=UNIFORM . F (IN DM IN (CUSTOMER)  , INDMAX (CUSTOMER) ,2) 
LET  SDELAY (COMPONENT)  =0 
FOR  1=1  TO  NPARTS  (COMPONENT)  ,  DO 
CREATE  A  REQN 


2QS)  “COMPONENT 
CSSDELAY 


LET  NHA  (EEC 
LET  IWAIT=I;  _  __ 

LET  WAIT  “IWAIT 

IF  IWAIT  LE  7  AND  IWAIT  GE  2.5 
SCHEDULE  AN  ISSUE  GIVEN  REON  IN  WAIT  DAYS 
ELSE  SCHEDULE  A  RECEIPT  GIVEN  REQN  IN  WAIT  DAIS 
REGARDLESS 
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LOOP 

SCHEDULE  AN  INDUCTION  GIVEN  CUSTOMER  IN  INDUCT/2  DAYS 
LOOP 

SCHEDULE  AN  END. SIM  IN  ENDSIM  DAYS 
SCHEDULE  A  START. STAT  IN  STARTSTAT  DAYS 
LET  K=KK 

START  SIMULATION 
LET  TIME. V=0 
REGARDLESS 
LOOP 
STOP 
END 

EVENT  DELIVERY 
LET  FLAG=0 
LET  TLOAD=N. QUEUE 
IF  K  EQ  9999 

SCHEDULE  A  DELIVERY  IN  N  DAYS 
REGARDLESS 

LET  TRANSCOST=TRANSCOST+S  HI PCOS  T 
FOR  EACH  REON  IN  THE  QUEUE,  DO 
REMOVE  THE  REQN  FROM  THE  QUEUE 
LET  COMPONENT=NHA  (REQN) 

LET  NP ARTS  (COMPONENTf=NP ARTS (  COMPONENT)  -1 
IF  NPARTS (COMPONENT)  LE  0 

LET  DLYlCOST (LINE  (COMPONENT))  =  (TIME. V-SDELAY (COMPONENT)  ) * 
COSTRATE  (LINE  (COMPONENT)  ) 

LET  NBRRFI  (LINE  (COMPONENT)  )=NBRRFI  (LINE  (COMPONENT))  +1 
REMOVE  THE  COMPONENT  FROM  REPAIR 
DESTROY  THE  COMPONENT 
REGARDLESS 
DESTROY  THE  REQN 
LOOP 
RETURN 
END 

EVENT  START. STAT 

FOR  EACH  CUSTOMER,  RESET  THE  TOTALS  OF  DLYlCOST 
FOR  EACH  CUSTOMER.  RESET  THE  TOTALS  OF  DLY2COST 
RESET  THE  TOTALS  OF  TLOAD 
LET  TRAN  SCOST=0 
FOR  EACH  CUSTOMER.  DO 
LET  NBRRFI  (CUSTOMER) =0 
LET  DLYlCOST  (CUSTOMER)  *0 
LET  DLY2COST (CUSTOMER) *0 
LOOP 
RETURN 
END 

EVENT  END.  SIM 

FOR  EACH  REQN  IN  THE  QUEUE,  DO 
REMOVE  THE  REQN  FROM  THE  QUEUE 
DESTROY  THE  REQN 
LOOP 

FOR  EACH  COMPONENT  IN  REPAIR.  DO 
REMOVE  THE  COMPONENT  FROM  REPAIR 
DESTROY  THE  COMPONENT 
LOOP 

LET  TC*0 

FOR  EACH  RECEIPT  IN  E V. S (I .  RECEIPT)  ,  DO 
CANCEL  THE  RECEIPT 
LOOP 

FOR  EACH  DELIVERY  IN  EV . S (I . DELIVERY)  ,  DO 
CANCEL  THE  DELIVERY 
LOOP 

FOR  EACH  ISSUE  IN  EV.  S  (I .  I S  SUE)  ,  DO 
CANCEL  THE  ISSUE 
LOOP 
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FOR  BACH  INDUCTION  IN  EV.  S  (I. INDUCTION)  ,  DO 
CANCEL  THE  INDUCTION 
LOOP 

SKIP  2  LINES 

PRINT  1  LINE  WITH  STARTSTAT,  ENDSIM  THUS 
STARTSTAT-***.  **  ENDSIM  =**  ** .  ** 

IF  K  EQ  9999 

PRINT  3  LINES  WITH  N,  TRANSCOST.  MN.N.QU,  AND  VAR.N.QU  THUS 
RESULTS  FOR  N=***.**  CASE.  TOTAL  TRANS  COST  IS  ****.** 

LOAD  MEAN=****.**  LOAD  VARIANCE****.** 
CUSTOMER  NBRRFI  AVE  LOC  D  COST  AVE  NON-LOC  D  COST 

ELSE 

PRINT  3  LINES  WITH  K ,  TRANSCOST,  MN.N.QU.  AND  VAR.N.QU  THUS 
RESULTS  FOR  K=**»  CASE.  TOTAL  TRANSPORTATION  COST*  ****.** 
LOAD  MEAN*****.**  LOAD  VARIANCE****.** 
CUSTOMER  NBRRFI  AVE  LOC  D  COSTS  AVE  NON-L  D  COSTS 

REGARDLESS 

FOR  EVERY  CUSTOMER.  DO 

PRINT  1  LINE  WITH  CUSTOMER.  N  BRRFI  (CUSTOMER)  , 

SDLY1COST  (CUSTOMER)  /(ENDSIM-STARTSTAT)  , 

SDLY2COST  (CUSTOMER) /(ENDSIM-STARTSTAT)  THUS 
***  ****  '****.*♦  ****.** 

LET  D 1 *SD LY1COST (CUSTOMER) +D1 
LET  D2=SDLY2COST (CUSTOMER) *D2 
LOOP 

LET  TC*  (TRANSCOST+DUD2)/  (ENDSIM-STARTSTAT) 

SKIP  t  LINE 

LET  PERsENDSIH-STARTSTAT 

PRINT  2  LINES  WITH  DI/PER.  D2/PER,  TC  THUS 

AVE  LOCAL  DEL  DELAY*  ****.**  AVE  OTHER  DELAY*  ****.** 

AVE  TOTAL*****. ** 

LET  TRANSCOST=0 

LET  D1 *0 

LET  D2=0 

RETURN 

END 


EVENT  INDUCTION  GIVEN  NCUST 

DEFINE  NCUST  AS  INTEGER  VARIABLE 

LET  CUSTOMER* NCUST 

CREATE  A  COMPONENT 

FILE  THE  COMPONENT  IN  REPAIR 

LET  LINE  (COMPONENT)  *CUSTOMER 

LET  Q= . 1 


LET  NPARTS  (COMPONENT)  -BINOMIAL. FfPMEAN  (CUSTOMER)  ,  0,4) 

LET  INDUCT=UNIFORH.  F (INDM IN?CUSTOMER)  ,INDMAX  (CUSTOMER)  ,2) 
LET  SDELAY  (COMPONENT! -TIME.  V  ♦  INDUCT/2 


LET  SDELAY  (COMPONENT) -TIME 
FOR  1=1  TO  NPARTS  (COMPONENT) 
CREATE  A  REQN 
LET  NHAjREQ 8)  -COMPONENT 
LET  IWA±T*lS5 DELAY 
LET  WAIT  -INDUCT/2+IWAIT 


IF  IWAIT  LE  7  AND  IWAIT  GE  2.5 

SCHEDULE  AN  ISSUE  GIVEN  REQN  IN  WAIT  DAYS 
ELSE  SCHEDULE  A  RECEIPT  GIVEN  REQN  IN  WAIT  DAYS 
REGARDLESS 
LOOP 

SCHEDULE  AN  INDUCTION  GIVEN  CUSTOMER  IN  INDUCT  DAYS 

RETURN 

END 


EVENT  ISSUE  GIVEN  ITEM 
DEFINE  ITEM  AS  INTEGER  VARIABLE 
LET  RBQN*ITEM 


FILE  THE  REQN  IN  THE  QUEUE 
IF  FLAG  EQ  0  AND  N. QUEUE  GE  K 
LET  FLAG*  1 

SCHEDULE  A  DELIVERY  NEXT 
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REGARDLESS 

LET  COMPONENT=NHA  (REQN) 

RETURN 

END 

EVENT  RECEIPT  GIVEN  DOC 
DEPINE  DOC  AS  INTEGER  VARIABLE 
LET  REQN=DOC 
LET  COMPONENT=NHA  (REQN) 

LET  NPARTS (COMPONENT) =NPART S (COMPONENT) -  1 
I?  NPARTS ( COMPONENT)  LE  0 

LET  DLY2COST  (LI  NE  (COMPONENT)  )  *  (TIME. V-SDELAY  (COMPONENT) )  *COSTRATE 
(LINE (COMPONENT)  ) 

LET  NBRHFI (LINE (COMPONENT) )=N BRRFI (LINE (COMPONENT)  )  *1 
REMOVE  THE  COMPONENT  FROM  REPAIR 
DESTROY  THE  COMPONENT 
REGARDLESS 
DESTROY  THE  REQN 
RETURN 
END 
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